A revisited history of cacao domestication in pre-Columbian times revealed by archaeogenomic approaches.

Humans have a long history of transporting and trading plants, contributing to the evolution of domesticated plants. Theobroma cacao originated in the Neotropics from South America. However, little is known about its domestication and use in these regions. In this study, ceramic residues from a large sample of pre-Columbian cultures from South and Central America were analyzed using archaeogenomic and biochemical approaches. Here we show, for the first time, the widespread use of cacao in South America out of its native Amazonian area of origin, extending back 5000 years, likely supported by cultural interactions between the Amazon and the Pacific coast. We observed that strong genetic mixing between geographically distant cacao populations occurred as early as the middle Holocene, in South America, driven by humans, favoring the adaptation of T. cacao to new environments. This complex history of cacao domestication is the basis of today's cacao tree populations and its knowledge can help us better manage their genetic resources.

Genome-wide association studies and genomic selection assays made in a large sample of cacao (Theobroma cacao L.) germplasm reveal significant marker-trait associations and good predictive value for improving yield potential.

A genome-wide association study (GWAS) was undertaken to unravel marker-trait associations (MTAs) between SNP markers and phenotypic traits. It involved a subset of 421 cacao accessions from the large and diverse collection conserved ex situ at the International Cocoa Genebank Trinidad. A Mixed Linear Model (MLM) in TASSEL was used for the GWAS and followed by confirmatory analyses using GAPIT FarmCPU. An average linkage disequilibrium (r2) of 0.10 at 5.2 Mb was found across several chromosomes. Seventeen significant (P ≤ 8.17 × 10-5 (-log10 (p) = 4.088)) MTAs of interest, including six that pertained to yield-related traits, were identified using TASSEL MLM. The latter accounted for 5 to 17% of the phenotypic variation expressed. The highly significant association (P ≤ 8.17 × 10-5) between seed length to width ratio and TcSNP 733 on chromosome 5 was verified with FarmCPU (P ≤ 1.12 × 10-8). Fourteen MTAs were common to both the TASSEL and FarmCPU models at P ≤ 0.003. The most significant yield-related MTAs involved seed number and seed length on chromosome 7 (P ≤ 1.15 × 10-14 and P ≤ 6.75 × 10-05, respectively) and seed number on chromosome 1 (P ≤ 2.38 × 10-05), based on the TASSEL MLM. It was noteworthy that seed length, seed length to width ratio and seed number were associated with markers at different loci, indicating their polygenic nature. Approximately 40 candidate genes that encode embryo and seed development, protein synthesis, carbohydrate transport and lipid biosynthesis and transport were identified in the flanking regions of the significantly associated SNPs and in linkage disequilibrium with them. A significant association of fruit surface anthocyanin intensity co-localised with MYB-related protein 308 on chromosome 4. Testing of a genomic selection approach revealed good predictive value (genomic estimated breeding values (GEBV)) for economic traits such as seed number (GEBV = 0.611), seed length (0.6199), seed width (0.5435), seed length to width ratio (0.5503), seed/cotyledon mass (0.6014) and ovule number (0.6325). The findings of this study could facilitate genomic selection and marker-assisted breeding of cacao thereby expediting improvement in the yield potential of cacao planting material.

Diversity and determinants of bitterness, astringency, and fat content in cultivated Nacional and native Amazonian cocoa accessions from Ecuador.

Cocoa (Theobroma cacao L.) is the only tree that can produce cocoa. Cocoa beans are highly sought after by chocolate makers to produce chocolate. Cocoa can be fine aromatic, characterized by floral and fruity notes, or it can be described as standard cocoa with a more pronounced cocoa aroma and bitterness. In this study, the genetic and biochemical determinants of sensorial notes and nonvolatile compounds related to bitterness, astringency, fat content, and protein content will be investigated in two populations: a cultivated modern Nacional population and a population of cocoa accessions collected recently in the Ecuadorian South Amazonia area of origin of the Nacional ancestral variety. For this purpose, a genome-wide association study (GWAS) was carried out on both populations, with results of biochemical compounds evaluated by near-infrared spectroscopy (NIRS) assays and with sensory evaluations. Twenty areas of associations were detected for sensorial data especially bitterness and astringency. Fifty-three areas of associations were detected linked to nonvolatile compounds. A total of 81 candidate genes could be identified in the areas of the association.

Population Structure of Moniliophthora perniciosa in the Main Cacao Producing Departments of Colombia.

The witches' broom (Moniliophthora perniciosa) is considered as one of the main threats for cacao production and, consequently, for chocolate production worldwide. In this work, the genetic diversity and population structure of M. perniciosa were analyzed for 59 isolates collected in five departments of Colombia and using 10 microsatellite markers. Analyses revealed 35 multilocus genotypes and clonal populations structure according to linkage disequilibrium analysis. One of the objectives of this study was to determine whether populations were differentiated by geographic origin or Theobroma cacao host genotype. Analysis of molecular variance, discriminant analysis of principal components, and Bruvo genetic distance suggested that the genetic structure was driven by geographic origin and not by T. cacao genotype. The results of this study were consistent with previous findings obtained in other cocoa-producing countries. Important insights were discussed regarding the dispersal patterns of the pathogen in Colombia and the genetic change of its populations because of different environmental conditions.

Integration of GWAS, metabolomics, and sensorial analyses to reveal novel metabolic pathways involved in cocoa fruity aroma GWAS of fruity aroma in Theobroma cacao.

Nacional is a variety of cocoa tree known for its "Arriba" aroma characterised mainly by fruity, floral, and spicy aromatic notes. In this study, the genetic basis of the fruity aroma of modern Nacional cocoa was investigated. GWAS studies have been conducted on biochemical and sensorial fruity traits and allowed to identify a large number of association zones. These areas are linked to both the volatile compounds known to provide fruity flavours and present in the beans before and after roasting, and to the fruity notes detected by sensorial analysis. Five main metabolic pathways were identified as involved in the fruity traits of the Nacional population: the protein degradation pathway, the sugar degradation pathway, the fatty acid degradation pathway, the monoterpene pathway, and the L-phenylalanine pathway. Candidate genes involved in the biosynthetic pathways of volatile compounds identified in association areas were detected for a large number of associations.

Two Main Biosynthesis Pathways Involved in the Synthesis of the Floral Aroma of the Nacional Cocoa Variety.

Theobroma cacao is the only source that allows the production of chocolate. It is of major economic importance for producing countries such as Ecuador, which is the third-largest cocoa producer in the world. Cocoa is classified into two groups: bulk cocoa and aromatic fine flavour cocoa. In contrast to bulk cocoa, fine flavour cocoa is characterised by fruity and floral notes. One of the characteristics of Nacional cocoa, the emblematic cocoa of Ecuador, is its aromatic ARRIBA flavour. This aroma is mainly composed of floral notes whose genetic and biochemical origin is not well-known. This research objective is to study the genetic and biochemical determinism of the floral aroma of modern Nacional cocoa variety from Ecuador. Genome-Wide Association Study (GWAS) was conducted on a population of 152 genotypes of cocoa trees belonging to the population variety of modern Nacional. Genome-Wide Association Study was conducted by combining SSR and SNP genotyping, assaying biochemical compounds (in roasted and unroasted beans), and sensory evaluations from various tastings. This analysis highlighted different areas of association for all types of traits. In a second step, a search for candidate genes in these association zones was undertaken, which made it possible to find genes potentially involved in the biosynthesis pathway of the biochemical compound identified in associations. Our results show that two biosynthesis pathways seem to be mainly related to the floral note of Nacional cocoa: the monoterpene biosynthesis pathway and the L-phenylalanine degradation pathway. As already suggested, the genetic background would therefore appear as largely explaining the floral note of cocoa.

Cacao breeding in Colombia, past, present and future.

Cacao (Theobroma cacao L.) is considered a key crop in Colombian social programs aiming at alleviating rural poverty, promoting peace in post-conflict regions and, replacing crops used for illicit purposes. Colombia is thought to be part of the center of origin of cacao; several germplasm collecting expeditions have been implemented, dating back to the 1940s. Despite that history, the first breeding program based on creating, selecting, and releasing full-sib progenies made extensive use of accessions introduced from other countries as parents. A new breeding strategy was adopted in the 1990s, based on mass selection of promising trees (high-yield and disease-resistant) in farmers' fields, resulting in the selection of clones released to farmers as planting material. In 2012, a new strategy, Recurrent Selection, was adopted by the Colombian Corporation for Agricultural Research, Agrosavia, based on the development of improved populations and allowing the selection of clones at the end of each cycle of recombination. The use of molecular markers is being integrated into this program in order to assist breeders in selecting material. This review provides details about the history and perspectives of the cacao breeding program in Colombia.

The use and domestication of Theobroma cacao during the mid-Holocene in the upper Amazon.

Cacao (Theobroma cacao L.) is an important economic crop, yet studies of its domestication history and early uses are limited. Traditionally, cacao is thought to have been first domesticated in Mesoamerica. However, genomic research shows that T. cacao's greatest diversity is in the upper Amazon region of northwest South America, pointing to this region as its centre of origin. Here, we report cacao use identified by three independent lines of archaeological evidence-cacao starch grains, absorbed theobromine residues and ancient DNA-dating from approximately 5,300 years ago recovered from the Santa Ana-La Florida (SALF) site in southeast Ecuador. To our knowledge, these findings constitute the earliest evidence of T. cacao use in the Americas and the first unequivocal archaeological example of its pre-Columbian use in South America. They also reveal the upper Amazon region as the oldest centre of cacao domestication yet identified.

Genomic preselection with genotyping-by-sequencing increases performance of commercial oil palm hybrid crosses.

BACKGROUND: There is great potential for the genetic improvement of oil palm yield. Traditional progeny tests allow accurate selection but limit the number of individuals evaluated. Genomic selection (GS) could overcome this constraint. We estimated the accuracy of GS prediction of seven oil yield components using A × B hybrid progeny tests with almost 500 crosses for training and 200 crosses for independent validation. Genotyping-by-sequencing (GBS) yielded +5000 single nucleotide polymorphisms (SNPs) on the parents of the crosses. The genomic best linear unbiased prediction method gave genomic predictions using the SNPs of the training and validation sets and the phenotypes of the training crosses. The practical impact was illustrated by quantifying the additional bunch production of the crosses selected in the validation experiment if genomic preselection had been applied in the parental populations before progeny tests. RESULTS: We found that prediction accuracies for cross values plateaued at 500 to 2000 SNPs, with high (0.73) or low (0.28) values depending on traits. Similar results were obtained when parental breeding values were predicted. GS was able to capture genetic differences within parental families, requiring at least 2000 SNPs with less than 5% missing data, imputed using pedigrees. Genomic preselection could have increased the selected hybrids bunch production by more than 10%. CONCLUSIONS: Finally, preselection for yield components using GBS is the first possible application of GS in oil palm. This will increase selection intensity, thus improving the performance of commercial hybrids. Further research is required to increase the benefits from GS, which should revolutionize oil palm breeding.

Deciphering the Theobroma cacao self-incompatibility system: from genomics to diagnostic markers for self-compatibility.

Cocoa self-compatibility is an important yield factor and has been described as being controlled by a late gameto-sporophytic system expressed only at the level of the embryo sac. It results in gametic non-fusion and involves several loci. In this work, we identified two loci, located on chromosomes 1 and 4 (CH1 and CH4), involved in cocoa self-incompatibility by two different processes. Both loci are responsible for gametic selection, but only one (the CH4 locus) is involved in the main fruit drop. The CH1 locus acts prior to the gamete fusion step and independently of the CH4 locus. Using fine-mapping and genome-wide association studies, we focused analyses on restricted regions and identified candidate genes. Some of them showed a differential expression between incompatible and compatible reactions. Immunolocalization experiments provided evidence of CH1 candidate genes expressed in ovule and style tissues. Highly polymorphic simple sequence repeat (SSR) diagnostic markers were designed in the CH4 region that had been identified by fine-mapping. They are characterized by a strong linkage disequilibrium with incompatibility alleles, thus allowing the development of efficient diagnostic markers predicting self-compatibility and fruit setting according to the presence of specific alleles or genotypes. SSR alleles specific to self-compatible Amelonado and Criollo varieties were also identified, thus allowing screening for self-compatible plants in cocoa populations.

Characterization of microsatellites identified by next-generation sequencing in the Neotropical tree Handroanthus billbergii (Bignoniaceae).

PREMISE OF THE STUDY: We developed microsatellite (simple sequence repeat [SSR]) markers in the Neotropical tree Handroanthus billbergii (Bignoniaceae), to be applied in assessment of genetic diversity in this species as a reference for inferring the impact of dry forest fragmentation in Ecuador. METHODS AND RESULTS: Using next-generation sequencing, we detected a total of 26,893 putative SSRs reported here. Using an ABI 3500xl sequencer, we identified and characterized a set of polymorphic markers in 23 individuals belonging to three populations of H. billbergii. CONCLUSIONS: We report a set of 30 useful SSR markers for H. billbergii and a large list of potential microsatellites for developing new markers for this or related species.

The coffee genome provides insight into the convergent evolution of caffeine biosynthesis.

Coffee is a valuable beverage crop due to its characteristic flavor, aroma, and the stimulating effects of caffeine. We generated a high-quality draft genome of the species Coffea canephora, which displays a conserved chromosomal gene order among asterid angiosperms. Although it shows no sign of the whole-genome triplication identified in Solanaceae species such as tomato, the genome includes several species-specific gene family expansions, among them N-methyltransferases (NMTs) involved in caffeine production, defense-related genes, and alkaloid and flavonoid enzymes involved in secondary compound synthesis. Comparative analyses of caffeine NMTs demonstrate that these genes expanded through sequential tandem duplications independently of genes from cacao and tea, suggesting that caffeine in eudicots is of polyphyletic origin.

The small RNA profile in latex from Hevea brasiliensis trees is affected by tapping panel dryness.

Natural rubber is harvested by tapping Hevea brasiliensis (Willd. ex A. Juss.) Müll. Arg. Harvesting stress can lead to tapping panel dryness (TPD). MicroRNAs (miRNAs) are induced by abiotic stress and regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs. This study set out to sequence miRNAs expressed in latex cells and to identify TPD-related putative targets. Deep sequencing of small RNAs was carried out on latex from trees affected by TPD using Solexa technology. The most abundant small RNA class size was 21 nucleotides for TPD trees compared with 24 nucleotides in healthy trees. By combining the LeARN pipeline, data from the Plant MicroRNA database and Hevea EST sequences, we identified 19 additional conserved and four putative species-specific miRNA families not found in previous studies on rubber. The relative transcript abundance of the Hbpre-MIR159b gene increased with TPD. This study revealed a small RNA-specific signature of TPD-affected trees. Both RNA degradation and a shift in miRNA biogenesis are suggested to explain the general decline in small RNAs and, particularly, in miRNAs.

The banana genome hub.

Banana is one of the world's favorite fruits and one of the most important crops for developing countries. The banana reference genome sequence (Musa acuminata) was recently released. Given the taxonomic position of Musa, the completed genomic sequence has particular comparative value to provide fresh insights about the evolution of the monocotyledons. The study of the banana genome has been enhanced by a number of tools and resources that allows harnessing its sequence. First, we set up essential tools such as a Community Annotation System, phylogenomics resources and metabolic pathways. Then, to support post-genomic efforts, we improved banana existing systems (e.g. web front end, query builder), we integrated available Musa data into generic systems (e.g. markers and genetic maps, synteny blocks), we have made interoperable with the banana hub, other existing systems containing Musa data (e.g. transcriptomics, rice reference genome, workflow manager) and finally, we generated new results from sequence analyses (e.g. SNP and polymorphism analysis). Several uses cases illustrate how the Banana Genome Hub can be used to study gene families. Overall, with this collaborative effort, we discuss the importance of the interoperability toward data integration between existing information systems. Database URL: http://banana-genome.cirad.fr/

BAC-end sequences analysis provides first insights into coffee (Coffea canephora P.) genome composition and evolution.

Coffee is one of the world's most important agricultural commodities. Coffee belongs to the Rubiaceae family in the euasterid I clade of dicotyledonous plants, to which the Solanaceae family also belongs. Two bacterial artificial chromosome (BAC) libraries of a homozygous doubled haploid plant of Coffea canephora were constructed using two enzymes, HindIII and BstYI. A total of 134,827 high quality BAC-end sequences (BESs) were generated from the 73,728 clones of the two libraries, and 131,412 BESs were conserved for further analysis after elimination of chloroplast and mitochondrial sequences. This corresponded to almost 13 % of the estimated size of the C. canephora genome. 6.7 % of BESs contained simple sequence repeats, the most abundant (47.8 %) being mononucleotide motifs. These sequences allow the development of numerous useful marker sites. Potential transposable elements (TEs) represented 11.9 % of the full length BESs. A difference was observed between the BstYI and HindIII libraries (14.9 vs. 8.8 %). Analysis of BESs against known coding sequences of TEs indicated that 11.9 % of the genome corresponded to known repeat sequences, like for other flowering plants. The number of genes in the coffee genome was estimated at 41,973 which is probably overestimated. Comparative genome mapping revealed that microsynteny was higher between coffee and grapevine than between coffee and tomato or Arabidopsis. BESs constitute valuable resources for the first genome wide survey of coffee and provide new insights into the composition and evolution of the coffee genome.

Identification of the Hevea brasiliensis AP2/ERF superfamily by RNA sequencing.

BACKGROUND: Rubber tree (Hevea brasiliensis) laticifers are the source of natural rubber. Rubber production depends on endogenous and exogenous ethylene (ethephon). AP2/ERF transcription factors, and especially Ethylene-Response Factors, play a crucial role in plant development and response to biotic and abiotic stresses. This study set out to sequence transcript expressed in various tissues using next-generation sequencing and to identify AP2/ERF superfamily in the rubber tree. RESULTS: The 454 sequencing technique was used to produce five tissue-type transcript libraries (leaf, bark, latex, embryogenic tissues and root). Reads from all libraries were pooled and reassembled to improve mRNA lengths and produce a global library. One hundred and seventy-three AP2/ERF contigs were identified by in silico analysis based on the amino acid sequence of the conserved AP2 domain from the global library. The 142 contigs with the full AP2 domain were classified into three main families (20 AP2 members, 115 ERF members divided into 11 groups, and 4 RAV members) and 3 soloist members. Fifty-nine AP2/ERF transcripts were found in latex. Alongside the microRNA172 already described in plants, eleven additional microRNAs were predicted to inhibit Hevea AP2/ERF transcripts. CONCLUSIONS: Hevea has a similar number of AP2/ERF genes to that of other dicot species. We adapted the alignment and classification methods to data from next-generation sequencing techniques to provide reliable information. We observed several specific features for the ERF family. Three HbSoloist members form a group in Hevea. Several AP2/ERF genes highly expressed in latex suggest they have a specific function in Hevea. The analysis of AP2/ERF transcripts in Hevea presented here provides the basis for studying the molecular regulation of latex production in response to abiotic stresses and latex cell differentiation.

Insight into the wild origin, migration and domestication history of the fine flavour Nacional Theobroma cacao L. variety from Ecuador.

Ecuador's economic history has been closely linked to Theobroma cacao L cultivation, and specifically to the native fine flavour Nacional cocoa variety. The original Nacional cocoa trees are presently in danger of extinction due to foreign germplasm introductions. In a previous work, a few non-introgressed Nacional types were identified as potential founders of the modern Ecuadorian cocoa population, but so far their origin could not be formally identified. In order to determine the putative centre of origin of Nacional and trace its domestication history, we used 80 simple sequence repeat (SSR) markers to analyse the relationships between these potential Nacional founders and 169 wild and cultivated cocoa accessions from South and Central America. The highest genetic similarity was observed between the Nacional pool and some wild genotypes from the southern Amazonian region of Ecuador, sampled along the Yacuambi, Nangaritza and Zamora rivers in Zamora Chinchipe province. This result was confirmed by a parentage analysis. Based on our results and on data about pre-Columbian civilization and Spanish colonization history of Ecuador, we determined, for the first time, the possible centre of origin and migration events of the Nacional variety from the Amazonian area until its arrival in the coastal provinces. As large unexplored forest areas still exist in the southern part of the Ecuadorian Amazonian region, our findings could provide clues as to where precious new genetic resources could be collected, and subsequently used to improve the flavour and disease resistance of modern Ecuadorian cocoa varieties.

Identification of novel microRNAs in Hevea brasiliensis and computational prediction of their targets.

BACKGROUND: Plants respond to external stimuli through fine regulation of gene expression partially ensured by small RNAs. Of these, microRNAs (miRNAs) play a crucial role. They negatively regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs (mRNAs). In Hevea brasiliensis, environmental and harvesting stresses are known to affect natural rubber production. This study set out to identify abiotic stress-related miRNAs in Hevea using next-generation sequencing and bioinformatic analysis. RESULTS: Deep sequencing of small RNAs was carried out on plantlets subjected to severe abiotic stress using the Solexa technique. By combining the LeARN pipeline, data from the Plant microRNA database (PMRD) and Hevea EST sequences, we identified 48 conserved miRNA families already characterized in other plant species, and 10 putatively novel miRNA families. The results showed the most abundant size for miRNAs to be 24 nucleotides, except for seven families. Several MIR genes produced both 20-22 nucleotides and 23-27 nucleotides. The two miRNA class sizes were detected for both conserved and putative novel miRNA families, suggesting their functional duality. The EST databases were scanned with conserved and novel miRNA sequences. MiRNA targets were computationally predicted and analysed. The predicted targets involved in "responses to stimuli" and to "antioxidant" and "transcription activities" are presented. CONCLUSIONS: Deep sequencing of small RNAs combined with transcriptomic data is a powerful tool for identifying conserved and novel miRNAs when the complete genome is not yet available. Our study provided additional information for evolutionary studies and revealed potentially specific regulation of the control of redox status in Hevea.

Discovery and mapping of a new expressed sequence tag-single nucleotide polymorphism and simple sequence repeat panel for large-scale genetic studies and breeding of Theobroma cacao L.

Theobroma cacao is an economically important tree of several tropical countries. Its genetic improvement is essential to provide protection against major diseases and improve chocolate quality. We discovered and mapped new expressed sequence tag-single nucleotide polymorphism (EST-SNP) and simple sequence repeat (SSR) markers and constructed a high-density genetic map. By screening 149 650 ESTs, 5246 SNPs were detected in silico, of which 1536 corresponded to genes with a putative function, while 851 had a clear polymorphic pattern across a collection of genetic resources. In addition, 409 new SSR markers were detected on the Criollo genome. Lastly, 681 new EST-SNPs and 163 new SSRs were added to the pre-existing 418 co-dominant markers to construct a large consensus genetic map. This high-density map and the set of new genetic markers identified in this study are a milestone in cocoa genomics and for marker-assisted breeding. The data are available at http://tropgenedb.cirad.fr.

Differential gene expression in different types of Hevea brasiliensis roots.

Three types of roots (taproots, first order laterals and second order laterals) were functionally characterized on 7-month-old in vitro plantlets regenerated by somatic embryogenesis in Hevea brasiliensis. A histological analysis revealed different levels of differentiation depending on root diameter. A primary structure was found in first and second order lateral roots, while taproots displayed a secondary structure. The expression of 48 genes linked to some of the regulatory pathways acting in roots was compared in leaves, stems and the different types of roots by real-time RT-PCR. Thirteen genes were differentially expressed in the different organs studied in plants grown under control conditions. Nine additional other genes were differentially regulated between organs under water deficit conditions. In addition, 10 genes were significantly regulated in response to water deficit, including 8 regulated mainly in lateral roots types. Our results suggest that the regulation of gene expression in lateral roots is different than that in taproots, which have a main role in nutrient uptake and transport, respectively.